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DECOUPLING IN UNEAR TIME -VARYING MULTIVARIABLE SYSTEMS 

By Viswanathan Sankaran* 

Langley Research Center 

SUMMARY 

This report deals with the problem of decoupling in linear multivariable systems 
which has received much attention recently. The necessary and sufficient condition for 
the decoupling of an m-input, m-output, linear time-varying dynamical system by state- 
variable feedback is described. The class of feedback matrices which decouple the sys- 
tem is characterized. Systems which do not satisfy these conditions are treated by par- 
titioning the system. Finally, necessary and sufficient conditions for decoupling a system 
which has both plant and observation noise are illustrated. 

INTRODUCTION 


Multivariable systems are defined as systems which have a multiplicity of inputs 
and outputs. Most of the problems associated with modern technology including aircraft 
and aerospace problems can be cast into such a framework. 


Multivariable systems are characterized by coupling or interaction which usually 
occurs between the variables of the system. Thus, in general, any given output is affected 
by variation of any of the input variables of the system. Similarly, variation of each input 
will affect all the output variables. Considerable effort has been directed recently toward 
restructuring the system so that the new system is noninteracting or decoupled. Since 
coupling is usually inherent in any plant, decoupling is a condition which is a part of the 
design objective. Thus, from an input-output point of view the new system possesses non- 
interaction, whereas from the point of view of original plant input and output, it is still 
interacting. 


Noninteracting multivariable systems can be conceived as consisting of a collection 
of individual subsystems, each of which has a single input and a single output. Dealing 
with single-input single-output subsystems has two advantages: (1) the problem of speci- 

fying performance requirements is simplified, and (2) each subsystem can be treated 
separately. With these advantages, the problem of designing to meet the needs is more 


niriUln Lnnmmn rt-f fL a win «tt noi fn/iliMi An rtf ottoi 1 o Vkl o ci o i ci ryrrl o rmf ni i f 

tFa.CLa.ux6, uctauDc ui tuc zaaoj uy it avunawic lv* uv^u». 


systems. 

*NRC-NASA Resident Research Associate. 
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Earlier work in formulating design procedures in which decoupling is required 
was reported by Voznesenskiy (ref. 1), Piven (ref. 2), Freeman (ref. 3), and Kavanagh 
(refs. 4 and 5). Horowitz (ref. 6) has discussed a practical method of synthesizing con- 
trollers to reduce the effect of plant parameter variations. Meerov (ref. 7) and 
Mesarovic (ref. 8) studied the structure of multivariable systems. All these contribu- 
tions use the transfer -function description of the dynamic systems. 

Morgan (ref. 9) in 1964 approached the problem based on the state variable repre- 
sentation. His main result was that the time invariant linear system 

x(t) = Ax(t) + Bu(t) 

y(t) = cx(t) 

can be decoupled if the matrix CB is nonsingular. Rekasius (ref. 10) extended Morgan’s 
result and outlined an essentially trial-and-error procedure for specifying a certain num- 
ber of system poles while the system is being decoupled. Falb and Wolovich (refs. 11 
and 12) in 1967 established the necessary and sufficient conditions for decoupling of 
multivariable systems by state variable feedback. They also made definite contributions 
to the synthesis problem. Gilbert (ref. 13) considered the complete structure of the solu- 
tion of the system. 

Porter (ref. 14), Viswanadham (ref. 15), Sankaran (ref, 16), and Majumdar and 
Choudhury (ref. 17) extended Falb and Wolovich's result to the time-varying case. 

Wonham and Morse (ref. 18) considered the problem in a geometric framework and they 
also proposed dynamic compensators (ref. 19) to achieve a desired pole distribution for 
the closed-loop-system transfer matrix. Extension of the problems with disturbances 
and problems with discrete time systems are considered by Sankaran and Srinath 
(refs. 16, 20 to 22). Decoupling using output feedback has been considered by Wang 
(ref. 23), Howze (ref. 24), and Singh and Rugh (ref. 25). The results for a class of non- 
linear systems have been reported by Nazar and Rekasius (ref. 26), Majumdar and 
Choudhury (ref. 27), and Tokumaru and Iwai (ref. 28). Application of decoupling theory 
to aircraft control problems can be seen in Cliff and Lutze (refs. 29 and 30) and in Yore 
(ref. 31). 

The purpose of this paper is to document the thesis work presented in reference 16 
which relates the requirements for decoupling multivariable systems by state variable 
feedback (refs. 11 and 12) to the time-varying case and includes a treatment of distur- 
bances in the systems. The emphasis is to give reasonable in-depth treatment of the 
modern problems facing the decoupling techniques and to present the results to date, but 
certainly not all the existing techniques for decoupling multivariable systems. The 
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results are based on Rozonoer's paper (ref. 32) on "A Variational Approach to the Prob- 
lem of Invariance of Automatic Control Systems.” However, this concept is extended to 
the multivariable decoupling problem. This report is organized in the following manner: 

First, the main results on the decoupling problem are presented. With the pre- 
liminary statement of the problem, the necessary and sufficient conditions for decoupling 
a linear dynamical system are given. The design approach for the feedback and feed- 
forward matrices is discussed and illustrated with an example problem. 

Secondly, systems which do not satisfy the given conditions for decoupling are 
treated and necessary and sufficient conditions to partially decouple these systems are 
given. An example is presented to clarify the results. 

Finally, systems with stochastic disturbances are considered. A concept of decou- 
pling is introduced and definitions and theorems applicable to such problems are 
presented. 

A bibliography of contemporary research works which concern decoupling is also 
given. In order to improve usability, the bibliography is arranged alphabetically. 

SYMBOLS 

A(t) system matrix of dimension n x n 

A c (t) compensated system matrix of dimension n x n 

B(t) input matrix of dimension n x m 

B c (t) compensated input matrix of dimension n x m 

Bj(t) jth column of matrix B(t) 

C(t) output matrix of dimension m x n 

C^(t) ith row of matrix C(t) 

D(t) decoupling matrix of dimension m x m 

D^t) ith row of matrix D(t) 

dj index parameter 
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E { } expected value operation on { } 

e i unit column vector with one in the ith element and zero elsewhere 

F(t) feedback matrix of dimension m x n 

G(t) feedforward matrix of dimension m x m 

H(t) matrix of dimension m x n 

I identity matrix 


m,n,sf 

indices 

J 4 (t) 

scalar function obtained as product C.(t) x(t) 

AJ(t) 

matrix with elements AJ. .(t) 

AJ U (t) 

increment in J.(t) due to any variation Acoj 

L(t) 

matrix of dimension m x n 

m j (t) 

arbitrary functions 

pHt) 

matrix of dimension n x m 

Q(t) 

matrix operator ^1 iL - A - BF^ 

Q k (t) 

operation of Q, k times 

s k,i (t) 

ith row of matrix S k (t) of dimension m x n 

s J,i (t) 

ith row of matrix S k c (t) of dimension m x n 

t,t*,r 

variable times 

to 

initial time 


(‘0 


= T = 


4 



final time 


u(t) m-vector representing the system input 

v(t) m-vector representing observation noise 

w(t) n-vector representing plant noise 

x(t) n-vector representing the state of the system 

x{t) minimum variance unbiased estimate of x(t) 

x(t) error in estimating x as x 

y(t) m-vector representing the system output 

Z. .(t) scalar function obtained as product d>.^( t) B. c (t) 

z(t) observation vector of order m 

ar^t) . . , cn^t) elements of P*(t) matrix 

<5( ) Dirac delta function 

A{t) arbitrary diagonal matrix of dimension m x m 

Xj(t) . . . A m (t) elements of A(t) matrix 

ix(t) scalar function which does not vanish simultaneously for any t 

0(t) adjoint system 

co(t) m-vector representing reference input 

A superscript T denotes the transpose of the quantity. A bar over a symbol used 
in partial decoupling represents partitioning of the quantity corresponding to the rank p 
of D(t). A tilda over a symbol in partial decoupling represents partitioning of the quan- 
tity corresponding to (m - p) terms. Dots over symbols denote derivatives with respect 
to time. The symbol V denotes "for all values." 
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PROBLEM STATEMENT AND RESULTS IN DECOUPLING 
BY STATE VARIABLE FEEDBACK 

Consider the continuous time linear dynamical system 


x(t) = A(t) x(t) + B(t) u(t) 

(1) 

y(t) = C(t) x(t) 

(2) 


where x(t) is a real n-vector representing the state of the system, and y(t) and u(t) 
are real m-vectors representing system output and input, respectively. The matrices 
A(t), B(t), and C(t) are n x n, n x m, and mxn matrices, respectively. Let F(t) 

and G(t) be mxn feedback and mxm feedforward matrices for all t e Jtg,t]j. 

A compensated system suitable for decoupling investigations using these feedback 
and feedforward matrices may be constructed by defining 

u(t) = F(t) x(t) + G(t) w(t) (3) 

The m-component cu vector is an external reference input to be used in controlling the 
m-component, y-output vector. (See fig. 1.) Different degrees of decoupling between the 
input and the output are defined later. 



Figure 1.- Block diagram for decoupled system. 

By using the state variable feedback, the compensated system can be written as 


x(t) = A c (t) x(t) + B c (t) co(t) 

(4) 

y(t) = C(t) x(t) 

(5) 
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where 


A C (t) = A(t) + B(t) F(t) 


B C (t) = B(t) G(t) 


It is assumed that the elements of A c (t), B c (t), and C(t) together with their first 
(n - 2), (n - 1), and (n - 2) derivatives, respectively, are continuous functions. 


Problem Statement 


Consider the function 


J.(t) = C i (t) x(t) 


(i = 1, 2, . . ,, m) (8) 


where C^(t) is the ith row of C(t). It follows that 

Jj(t) = y ± (t) (' 

Let AJ. .(t) be the increment in «L(t) due to any variation Aon in cu.(t) 
i,J 1 J J 

(t Q = t = t). Also let AJ represent a matrix with elements AJy. 

Definition 1: The system given by equations (4) and (5) is said to be uniformly 

decoupled for all t e [to,ti] if (1) the variation in y.(t) due to any variation in ^(t) 

^ £ r s t) whenever i # j is identically zero V t e and (2) the variation in 

y^(t) due to any variation in = T = tj is not zero V te From this 

definition, it follows that 


AJ..(t) = 0 


AJ t .(t) * 0 


V t e t 0 ,t]]; i * jj (10) 


V t e [tg^j ; i,j = 1 , 2 , . . mj (11) 


for uniformly decoupled systems. 
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Main Results for Decoupling by State Variable Feedback 

Theorem 1 : For AJ^ j(t) to be identically zero on jtQ,tij, 
i,j = 1, 2, . . m; i * j, it is necessary and sufficient that 


C i (t) Q k (t) B°(t) = 0 


h ^[Wi} i ^ 

ij = 1 , 2, . * m; 

U = 0,1, . . ., n-1 J 


( 12 ) 


where the matrix operator Q is given by 


Q = l|-A c 


(13) 


Q k denotes the operation of Q k times and B. c (t) is the jth col- 
umn of B c (t) (ref. 32). J 

Proof: The quantity AJ. ^(t), the variation in J. (i = 1, 2, . . m) due to a vari- 

ation in Wj(r) = t = t; i * j) may be shown to be given by 


AJ i,j(t) = - 0 i T ('!’) B. c (t) AcDj(t) dr 

where $.(r) satisfies the adjoint equations 
^ r) = -A cT (t) 0.(t) 

and 

^i(t) = -C. T (t) (i = 1, 2, . . ., m) 


(14) 


(15) 


(16) 


where the superscript T denotes the transpose of the quantity. 
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Necessity: From equation (14), it is stated that equation (10) is satisfied if and 
only if 


Z y (r) = 0. T (r) Bj C (t) = 0 



If equation (17) is fulfilled, then the fulfillment of equation (10) is clear by virtue of 
equation (14). The converse is true, namely, equation (10) implies equation (17). To 
prove this statement, assume that equation (17) is not fulfilled at some single point 

t e |tQ,tJ, for example suppose that </).^(t) B. c (t) > 0. Then the continuity of the func- 
tion Z. .(t) implies the existence of an entire segment including the point t on which 

1 ?J 

T 

<p . (r) B, c (t) > 0. Now by choosing the increment Aw.(r) equal to zero outside this 

segment and equal, for example, to unity inside it, it can be seen that by virtue of equa- 
tion (14), AJ. .(t) < 0, that is, equation (10) is not satisfied which is a contradiction. In a 
1 ?J 

similar argument, it can be shown that equation (10) is not satisfied for <^k (t) Bj (r) < 0. 
Thus the statement has been proved. 

Differentiating equation (17) and making use of equation (15) and the definition for 
Q result in 


d K Z. .(r) 


±f— = V M QV) B. C (t) = 0 

k * ■> 


dr 


i * j; i,j = 1, 2, . 
\k = 0, 1, . . ., n-1 


(18) 


Substitution of equation (16) into equation (18) yields 




dr* 


= -C.(t)[Q k (r) B. c (t)1 = 0 

T~t L JT=t 


A * ), i; j - 1, 2, . . ., m) 

\k = 0, 1, . . n-1 / 


(19) 


Since t is arbitrary, equation (12) follows and hence the necessity is proved. 

Sufficiency: From equation (12), because C (t) is not identically zero, it follows 
that the vectors Bj C (t), QBj C (t), . . ., Q n-1 Bj C (t) are linearly dependent; that is, there 
exist functions ju^t), pij(t), . . ., which do not vanish simultaneously for any 

t for which 


n-1 

£ Mg(t) Q S Bj C (t) = 0 
s=0 


( 20 ) 
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From equations (18) and (20) it follows that the function -(t) satisfies the differential 
equation 


n-1 

s=0 at 


0 


with the boundary condition given by equation (19). The unique function which satisfies 
this relation is 


Zf j(t) = tpjit) B. c (t) = 0 


and hence the sufficiency follows. 



To apply the results of theorem 1 to the synthesis problem, it is convenient to 
express equation (12) in terms of the row vectors s£ ^(t) (k = 0, 1, . . n - 1; 
i = 1, 2, . . ., m) defined by 




So4<‘) - c,(t) 

&14 (t) =S s k4 (t, + ^ <t)AC(t) 




( 21 ) 


Lemma 1 : Theorem 1 holds if equation (12) in the statement of the 
theorem is replaced by 


S k4 (t > B j C « = 0 


Proof: From equation (21), 


(v ‘ « [Vi]i » J:\ 


ji,j — 1,2,. . . , m; 
\ k = 0, 1, . . ., n-1 


I 


( 22 ) 


C.(t) Q k (t) B. C (t) = Y (-l) m i^Fs^ ( t ) B . c (t) 

^ J , , K" 1 m niji j 


m=0 



( 23 ) 


From the results of the theorem 1 and equation (23), the lemma is proved. 
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To obtain the inequality condition of equation (11), the concepts of uniform output 
controllability or uniform input observability are used. By use of the dual relation pre- 
sented by Silverman and Meadows (ref. 33), the following comments can be made: 

Comment 1: The ith subsystem is uniformly input observable on jto,tjJ if and only 
if the rank of the vector 






(24) 


is one for all 



From the definition of uniform decoupling, if the system given by equations (4) 
and (5) satisfies simultaneously equation (22) and comment 1, then it is uniformly 
decoupled. 

From the preceding discussion and theorem 1, it follows that if the n x m matrix 


P\t) = 


S C Q>i (t) B c (t) 
Sj 4 (t) B c (t) 


(i = 1, 2, . . ., m) 


(25) 


s'.i ;i (t) B c (t) 


i ! 

has all the columns identically zero except the ith column, that is, it has at least one non- 
zero element 


P*(t) = 


^(t) 

«o l (t) 


0 0 


0 0 


(i = 1, 2, . . m) 


“nW 


then it is uniformly decoupled. 


(26) 
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Comment 2: The system is uniformly decoupled for all t e [to,tiJ, if on ty if 
for all i = 1, 2, . . m, P*(t) is of the form of equation (26) so that the rank of its ith 
column is 1 for all t e |tg,t]_J. t 


Feedback and Feedforward Matrices 


In order to obtain the requisite F(t) and G(t) consider the indices dj defined 
as 

dj A Min|j: S. .(t) B(t) * 0; V t e [tg,tjj; j = 1, 2, . . n-lj (27) 

where the row vectors S, .(t) (j = 0, 1, . . n - 1; i = 1, 2, . . m) are defined by 

J J 1 

S 0,i (t > ■ c i“) 1 

) (28) 

W B -SV‘ , *V» AW 

J 

The dj values are, in general, time dependent. But if S. .(t) B(t) ± 0 for some 

* J ’ 1 

j and some t = t , then the continuity and differentiability assumptions on A(t), B(t), 

and C(t) imply that Sj ^(t) B(t) * 0 on some subinterval including the point t*. Hence, 
when dj values exist, they are constants over some subinterval of £tg,t^J. In the mate- 
rial that follows, it is assumed that the dj values are constants over [to,tiJ. 

Lemma 2 ; s£ ^(t) = ^(t) for k = dj and for all t € jtQ,t^. 


Proof: From equation (21) and the definition of Sq .(t) 

syt) - c,® = s 01 <t) 

the proof can be completed by induction. Assume for some p < d^ 


Sdpp.lM - W° 


V t € 


[Wx] 
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Then for - p + 1 by equation (21), 


s 5 i - P+ i,i (t) =|[ s 5 i -p.i (t> ] + s 5 i -p4 lt)AC(,) 


From equations (30) and (6), 


-^W*] + V-‘ (t) A<t> + Sd i-P' 1<l) B(t> F<t> 


S V P+ l,i (t) ' S d, W‘> + V/ 0 B(,) F(t) 


In view of equation (27) and the fact that d^ - p < dj, 


£ Ve*i,i<‘> - W 


Thus, by induction it follows that 


S M«> * V*’ 


V t £ r Vl ]; k S dj 


Hence the lemma is proved. 


Di(t) = S d ..(t) B(t) 


(i = 1, 2, . . m) 


and let D(t) be the matrix whose ith row is given by equation (36); that is, 


s di,i<‘> B <‘> 


s d 2 ,2 (t) B(t) 


D(t) = 


S d (t) B(t) 

u m> m 

u - 
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Let A(t) be an arbitrary diagonal matrix, nonsingular for all t e [to,ti"L namely, 


A x (t) 


x 2 (t) 


x 3 (t) 


A(t) = 


( 38 ) 




and let 


S d 2 +l,2 (t > 


H(t) = 


5 d m+ l,m (t) 


(39) 




Theorem 2 : The system given by equations (4) and (5) is uniformly 
decoupled if and only if D(t) is nonsingular for all t e jtQ.tjJ: Then 

the class of decoupling feedback matrices is given by 


F(t) = -D -1 (t) H(t) 

(40) 

G(t) = D -1 (t) A(t) 

(41) 


Proof: 

Sufficiency: In view of equations (27) and (35), equation (25) can be written as 
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0 


PV = 


sS .(t) B (t) G(t) 
u i ?i 


S n-l,i (t) B(t) G(t) 


(i = 1, 2, . . m) (42) 


Now 


S^.(t) B(t) G(t) = S di?i (t) B(t) G(t) = D.(t) D _1 (t) A(t) = e. T A(t) (v t e (43) 


S d i ,i (t) B(t) G(t) = Vi 


V t e 




(44) 


where is a unit row vector with one in the ith place and zero elsewhere. From 

equations (36) and (40) 


s d i+ l,i (t) ■ Vu (t) + V * 1 B(t) F(t) ■ - D i (t) D " 1(t) H(t> - 0 

(v t € to.tjJ; i = 1, 2 raj (45) 


From equations (45) and (21), it follows that 


s W t} = 0 


V t e 


[Vi]; 


k S 1 
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Hence, 


pi(t) = 


0 0 A.(t) 0 ... 0 


(46) 


Since A^(t) 4- 0 for all t e P*(t) has the form of equation (26) with the rank of 

its ith column 1 for all t. Hence by comment 2 the system is uniformly decoupled. 

Necessity: Suppose the system is uniformly decoupled, then by comment 2, pi(t) 
is of the form of equation (26), so that for each t e jtQjtjl, at least one member of its 
ith column is nonzero. From equations (27) and (41), it follows that D^L) G(t) = A^t) e. T 
where A^t) * 0 for all t e Jt^tjJ (i = 1, . . ., m). 

Thus, 


D(t) G(t) = diag ^A^(t),A, 


(t),. 






A 


(47) 


The matrix G(t) is desired to be nonsingular for all t e jjtQ,tjJ in order that all 

the inputs influence all the outputs. Hence it follows from equation (47) that D(t) is non- 
singular and hence the theorem. 

The matrix F(t) obtained in equation (40) is a class of matrix to insure that P*(t) 

has rank one for i = 1, 2, . . ., m and for all t e |tQ,tjJ. As is evident, there is no 

choice in determining F(t). However, by assigning arbitrary coefficients, a particular 
form of F(t) still insuring that the rank of PMt) is one is given by the following 
comment: 
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Comment 3: The rank of P*(t) is one for all t € 
F(t) = -D -1 (t) L(t), where 


A 

S d 1 + M (t) + X m j 1(t) S j,l (t) 

j=0 


L(t) = 


m 


- •* 

t 0 ,t x (i = 


= 1, . . m) if 


W* + I m ! m<t) 

j=o 


(48) 


where the m^(t) coefficients are arbitrary coefficients available to the analyst to insure 
that the compensated system is stable. 

As is evident H(t) of equation (39) is a particular case of equation (48). 

An Example Problem 

Consider the system given by equations (1) and (2) where 


A(t) = 


1 4 + cos t 

4 + cos t sin t 


0 

2 


(49) 


B(t) = 


e- 1 

1 

0 


0 

0 

4 + cos t 


1 

0 

1 


(50) 
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and 


C(t) = 


4 + cos t 
0 
0 


0 

0 

4 + cos t 


(51) 


Then from equation (37) 


Now 


D(t) = 


S 0>1 (t) B(t) 
S 0 2 (t) B(t) 
S 0,3 ( t) 


D(t) = 


(4 + cos t)e _t 0 

0 4 + cos t 

4 + cos t 4 + cos t 


D _1 (t) = 



1 

o 

1 

-1 

4 + cos t 

4 + cos t 


1 


4 + cos t 

1 

1 


-1 


,-t 


4 + cos t 
e - t 


1 

.-t 


4 + cos t 
-e _t 


(52) 


(53) 


Because D is nonsingular, the system can be decoupled and the corresponding feedback 
and feedforward matrices are given by equations (40) and (41); that is, 


F(t) = 


-1 


4 + cos t 


11 

f 12 

f 13 

21 

f 22 

f 23 

31 

f 32 

f 33 


(54) 
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and 


where 


G(t) = 


-X 


1 


4 + cos t 4 + cos t 


-X 2 


, X 2 e 

Art — ' “ 

“ 4 + cos t 


x 2 e 


-t 


X 3 e 


-t 


4 + cos t[ 

-t 


-X 3 e 


(55) 


f 1 1 = (4 + cos t)* 


f 12 - (4 + cos t)(sin t + m n ) - sin t 


f - A j. m ^ ™ 2 

f 13 ~ 4 + 0 " 0 


f 01 = -(1 + e -t + mA) + — 1 + e *(4 + cos t) 

' u / 4 + cos t 


. lm n + 2 sin ti / 3 

f 00 = -e' 1 - (4 + cos t) + \sin t + idq 


22 


4 + cos t 


m n ^e‘ t mAe -t 
f = -i + — y — + — y 

4 + cos t 4 + cos t 


f 3 1 = -sin t + (4 + cos t) 


^1 + e -t + mQ 1 


+ (4 + cos t)e 


-t 


f 32 " e 


-t 


1 / 3 

niQ + 2 sin t - (4 + cos t)fsin t + niQ 


+ (4 + cos t)' 


t/ 2 3 

f 33 = -(4 + cos t) + e 'I -4 + m Q - m Q 
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and the m. 1 and Xj terms are arbitrary and may be used to insure stability of the 

compensated system. The compensated matrices can be obtained by using equations (6) 
and (7) and these relations; that is, 


B c (t) = 


4 + cos t 




4 + cos t 4 + cos t 


and 


(56) 


A c (t) = A(t) + B(t) F(t) 

Time -Invariant System 

The preceding results are available for time -invariant systems (ref. 12). The nec- 
essary and sufficient condition for decoupling by state variable feedback is that the 
matrix D whose ith row is given by 

di 

= C^A B (57) 


C^B * 0; j = 0, 1, . . ., n-lj = n - 1 

(c^B = 0 for all jj (58) 

The decoupling feedback and feedforward matrices F and G are given by 

F = -D -1 A* (59) 


be nonsingular, where 


dj = Min^j: 


and 


G = D 


-1 


(60) 
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where 



( 61 ) 


PARTIAL DECOUPLING BY STATE VARIABLE FEEDBACK 

In the previous sections it was shown that for decoupling by state variable feedback, 
it is necessary and sufficient that the decoupling matrix D(t) is nonsingular. When such 
is not the case, it is desirable to decouple at least part of the output from part of the ref- 
erence input. 

Problem Statement 

Consider for the system given by equations (1) and (2). Let the rank of 
D(t) = p < m. Partition the output vector as 

y 

y = 

y 


where y is a p-vector corresponding to the rank of D(t) and y, a m - p vector. 
Similarly, partition the input vector as 

u 

u = — 
u 

where u is a p-vector and u, a m - p vector. 
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Then equations (1) and (2) can be written as follows: 


k(t) = A(t) x(t) + [B(t);B(t)] 


u(t) 

u(t) 


( 62 ) 


y(t) 


! oi 
1 <-+ 

' w ■ 

y(t) 


c(t) 


x(t) 


(63) 


f+wt mm*— ~ 

where B, B, C, and C are matrices of compatible order. Let the reference input 
vector be w corresponding to the p inputs u. Let F(t) and G(t) be pxn and 
p x p matrices for all t e |tg,t jJ. Then let the input u to the system using state vari- 
able feedback be (fig. 2) 


G(t) = F(t) x(t) + G(t) a>(t) 


(64) 



Figure 2.- Block diagram for partially decoupled system. 

Consider the function 

j.(t) =C t (t) x(t) (i-1,2, . . ,,p) (65) 

J 4 (t) = C.(t) x(t) (i = 1, 2, . . ., m-p) (66) 

where Cj(t) is the ith row of C(t) and C^(t) is the ith row of C(t). From equa- 
tion (52) it follows that J^t) = y^t) and Jj(t) = y^t). 
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Let AJ. .(t) be the increment in J,(t) due to any variation Awj in u>.(r) 

1)J 1 J J 

= t = tj and AJ^ j(t) be the increment in J^t) due to any variation AWj in Wj(r) 

= t = t^. Also let AJ and AJ represent matrices with elements AJ i j and AJ^ j, 
respectively. 

Definition 2: The system given by equations (62) and (63) is said to be uniformly 
partially decoupled for all t e ftQ,^ if ( 1 ) the variation in y^t) due to any variation 
in ^j( T ) ^0 = r = t; j * ij is identically zero V t e Jt 0 ,t J, (2) the variation in y^(t) 

due to any variation in aL(r) |tg = t I tj is not zero V t e an d (3) ||AJ(t)|| is 


minimized. 

From this definition, it follows that 


Hj(t) = o 




(67) 



and 



( 68 ) 



m-p 

- s "p 1 

* £ M 1 



4, u (t > 


(69) 


is minimized. 


Theorem 3 : For the system given by equations (67) and ( 68 ) to be uni- 

formly partially decoupled in the sense of the preceding definitions, 

(1) It is necessary and sufficient that 


LjW vit] B.(t; = u 


(y t C [Vi]; 

i,j = i, 2 , . . p; I 
\k - 0 . 1 , . . ., n -1 j 


(70) 
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where Q is the matrix operator |l - A - BfJ and 

(2) To minimize the cross -coupling effect, it is sufficient that 


m-p p 


sup y y c^t) B.w 

M iri i=1 


(k = 0, 1, . . n-1) (71) 


is minimized. 


Proof: The proof follows along the lines of theorem 1 and can be found in refer- 
ence 16. 

Comment 4: Equivalently, for the system to be uniformly partially decoupled, 

(1) it is necessary and sufficient that D(t) whose ith row 

^(t) = S d . 4 (t) B(t) (i = 1, 2, . . ., p) (72) 

is nonsingular |V t € fto,tj"]^ and (2) to minimize the cross-coupling effect it is sufficient 
that 

D(t) 


S d .(t) B(t) 


(i = p + 1, 


m) 


(73) 


is minimum where the norm is selected as given by equation (69). 

In order to relate these results to the D(t) of the original system (1) and (2), par- 
tition D(t) as follows: 


D(t) ! D(t) 


D(t) 


I 

-*V 


(74) 


D(t) ! D(t) 


where D(t) is a pxp matrix of rank p t e |tg,t]j\ D(t) is ra-pxm-p 
*“ / 
matrix, D(t) is m - pxp matrix, and D(t) is p x m - p matrix. The matrix D(t) 
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is invertible and hence the design procedure for F(t) and G(t) is identical to ordinary 
decoupling. 

It may be noted that the rows of D(t) correspond to the output y(t) and the col- 
umns of D(t) correspond to the input u(t). This fact gives us an idea of what outputs 
and inputs to choose for decoupling purposes. 

Comment 5: The classes of decoupling matrices of equation (64) are given by 


F(t) = -D -1 (t) L(t) 

(75) 

G(t) = D _1 (t) A(t) 

(76) 


where 


L(t) = 


l l 

s d, + l,x (t) * 1 “j 1(t) s i,l (t) 

1 1=0 


p j=0 


(77) 


where 



terms are arbitrary coefficients and A(t) 


trary coefficients and is nonsingular 



is a diagonal matrix with arbi- 


An Example Problem 

As an example, consider the system given by 


A(t) = 


4 + cos t 

-1 


4 + cos t 
sin t 
0 


-1 

0 

2 


(78) 
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0 


1 


e " 1 

B(t) =0 0 0 

0 4 + cos t 1 


( 79 ) 


4 + cos t 


C(t) = 


0 

0 



0 

4 + cos t 


0 

1 

1 


Then 


(4 + cos t)e 


-t 


D(t) = 


0 

0 


0 

4 + cos t 
4 + cos t 


4 + cos t 
1 
1 


(80) 


(81) 


and D(t) is singular for all t e [0,°°). Hence the system cannot be decoupled. How- 
ever, for t e Co,°°), the rank of D(t) is two and so two of the outputs can be decoupled 
from the reference inputs and the cross- coupling effects minimized by using these results. 

Without loss of generality, consider the outputs to be decoupled to be y^ and y2* 

Then, 


C(t) = 


4 + cos t e -t 0 


0 1 


and 


C(t) = [o 4 + cos t lj 


(82) 


In order to minimize D in the sense of equation (69), that is, the cross coupling between 
the output y 3 and the reference inputs and assign the inputs u^ and u 3 to 
and WgJ that is, 


26 



(83) 


(84) 


(85) 


( 86 ) 


(87) 


( 88 ) 


(89) 


27 



where 


G(t) = ( 4 + cos ^ e_t 6 * 


f t1 = < (4 + cos t) 2 + e~* + m n ^(t) - sin t>e 


fl2 = + sin t - 1 + e*(4 + cos t) 2 


fig = -e l (4 + cos t) 3 + m0 2 (t) 


f 91 = -(4 + cos t) 


f 22 " 0 


■^23 = ^ + m Q (t) (4 + cos t) 


and the m. 1 and X. terms are arbitrary and may be used to insure stability of the par- 

J J 

tially compensated system. The corresponding compensated system matrices are given 
by 


a* i a. in 0 


11 

a 12 

c 


21 

a 22 

0 

0 


28 



where 




sin t 
4 + cos t 




1 • * 

+ sin t - 

4 + cos t 


i) 


-t 

e 



= 4 + cos t 



= sin t 



and 


0 


* 2 <t)J 

STOCHASTIC DECOUPLING BY STATE ESTIMATOR VARIABLE FEEDBACK 

In the previous sections decoupling of linear multivariable system by state feedback 
and a transformation of the input is considered. Even though output feedback can be used 
for decoupling, some of the flexibility obtained with state variable feedback as in specify- 
ing closed loop poles, and so forth, will, in general, be lost. If the states are not directly 
accessible, a possible approach is to construct an observer to reconstruct the inaccessible 
states to generate the control law. In this section stochastic decoupling in linear multi- 
variable system with known plant and observation noise is considered. A definition of 
decoupling applicable to this and a procedure for achieving such decoupling are presented. 


B c (t) 


*l<t) 


4 + cos t 
0 
0 


(94) 


Consider the linear plant and observation model 
x(t) = A(t) x(t) + B(t) u(t) + w(t) 
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y(t) = c(t) x(t) 


(95) 


z(t) = y(t) + v(t) (96) 

where x(t) is a real n-vector representing the state of the plant, y(t), u(t), and z(t) 
are real m-vectors representing the plant output and input and observation, respectively. 
A(t), B(t), and C(t) are matrices of compatible dimensions. The vectors w(t) and v(t) 
are plant and observation noise vectors that are assumed to be zero mean and Gaussian 
with 

Cov{w(t),w(r)} = Q(t) 6(t - r) 

Cov{v(t),v(r)} = R(t) 6(t - r) > (97) 

Cov{w(t),v(r)} = S(t) 6(t - r) 

and x(tg) is a Gaussian random variable with known mean and covariance which is inde- 
pendent of w(t) and v(t) for all t = tg. Let x(t) be the minimum variance estimate 
of x(t), and let the control law be given by 

u(t) = F(t) x(t) + G(t) w(t) (98) 


where the m-vector o>(t) is the reference input and F(t) and G(t) are matrices of 
compatible dimensions. 

Consider the random functions, 

J.(t) = C^t) x(t) (i = 1, 2, . . ., m) (99) 


where Cj 
variation 



is the ith row of C(t). 
Acuj in u^(r) ^tg = t = 


Let AJ. >;j (t) 
t) and let’ AJ 


be the increment in J^t) due to any 
represent the matrix with elements 


Definition 3: The system given by equations (94) to (96) is said to be stochastically 
decoupled, if (1) the expected variation in y.(t) due to any variation in w.(r) ^ S t ^ tj 

if i * j is zero; and (2) the expected variation in y.(t) due to any variation in co. (r) 

(tg = r — t) is not zero. 1 1 


Mathematically, 


E 



= 0 


(i * j) (100) 
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and 




. ., m; V t € 



( 101 ) 


Theorem 4 : For the stochastic control problem to satisfy the equality 
relation of equation (88), it is necessary and sufficient that 


C t (t) Q k (t) Bj(t) = 0 


* j; i,J = 1, 2, . . 

k = 0, 1, . . n-1; 



( 102 ) 


where 


Q 


is the matrix operator 




and 


B i<t) 


is the jth 


column of B(t). 


Proof: Because 


x = x - x (103) 

where x is the error in estimating x as x, then 

x(t) = (A + BF)x + BGco + w - BF2 (104) 

If x represents an unbiased estimate, 

E{x(t)} = 0 

The rest of the proof follows along the lines of theorem 1. 

Comment 6: For the stochastic control problem to be decoupled, it is necessary 
and sufficient that D(t) whose ith row 

Dj(t) = S^it) B(t) (i = 1, 2, . . m) (105) 

is nonsingular t e 



31 



The classes of feedback and feedforward matrices are the same as for ordinary 
decoupling. Also the results can be extended in a straightforward manner to the case 
when the D matrix is singular to obtain a partially decoupled system. 

CONCLUDING REMARKS 

The necessary and sufficient conditions for decoupling, time-varying, continuous- 
time, linear-multivariable systems have been presented in this paper. A constructive 
method for determining the required feedback and feedforward matrices was described 
and illustrated in an example problem. 

For systems which do not satisfy the necessary and sufficient conditions, a method 
for accomplishing partial decoupling was presented and illustrated through an example 
problem. Finally, a definition of decoupling applicable to linear systems with known plant 
and observation noise and a procedure for decoupling such systems were given. 

Although the procedures discussed in this report apply to linear continuous-time 
systems, it is obvious that with minor modification they would apply equally well to 
discrete-time systems. 

In spite of the generality of the present theory, there remain a number of important 
problems for future investigation. Probably the most challenging of these is to extend the 
present theory to a broader class of systems. Although efforts have been made to develop 
a decoupling theory for nonlinear systems and distributed systems, the considerable dif- 
ficulties of these generalizations of the problem so far have inhibited substantial progress. 

Another important endeavor would be to develop meaningful questions regarding the 
sensitivity of decoupled systems. Efforts in this direction were made by the author and 
presented in the Proceedings IEEE Trails., Feb. 1973 (pp. 241-242). In fact, it is clear 
that in practice no physical system can be decoupled exactly. Obviously, what is needed 
is a useful and precise definition for an ''approximately” decoupled system. Quantitative 
considerations such as this one are of considerable practical importance and almost 
nothing is known about them. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., October 12, 1973. 
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